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The sensitivity of elastic breakup to the interior of the projectile wave function is analyzed. 
Breakup calculations of loosely bound nuclei (*B and '^^Be) are performed with two different de- 
scriptions of the projectile. The descriptions differ strongly in the interior of the wave function, but 
exhibit identical asymptotic properties, namely the same asymptotic normalisation coefficient, and 
phase shifts. Breakup calculations are performed at intermediate energies (40-70 MeV/nucleon) on 
lead and carbon targets as well as at low energy (26 MeV) on a nickel target. No dependence on the 
projectile description is observed. This result confirms that breakup reactions are peripheral in the 
^— ^ ■ sense that they probe only the external part of the wave function. These measurements are thus 

' not directly sensitive to the total normalization of the wave function, i.e. spectroscopic factor. 
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I. INTRODUCTION 

(N I 

^ I The development of radioactive ion beams in the mid-80s has allowed the study of nuclei far from stability. This 
technical breakthrough led to the discovery of halo nuclei on the neutron-rich side of the valley of stability pj, i^. 
These loosely bound nuclei have a strongly clusterized structure [3, 0, B @] ■ In a simple model, they are seen as a 
core, that contains most of the nucleons, to which one or two neutrons are loosely bound. Due to this small binding, 
the valence neutrons tunnel far outside the classically allowed region and form a sort of halo around the core [3]. 
\^ , Although less probable, proton halos are also possible. 

• Being very short lived, these nuclei can be studied only by indirect techniques like elastic breakup. In this reaction, 
' the halo dissociates from the core through interaction with a target [1, 0, H, ■ From the comparison of the data 
> • with the theoretical prediction, one usually extracts a spectroscopic factor for the single-particle configuration of the 
1— I [ projectile model [1, B [S El • Alternatively, it has been suggested that breakup reactions being highly peripheral, 
^ . only the asymptotic normalization coefficient (ANC) can be extracted from the measurements fl2l. Iisj. These two 
viewpoints are contradictory and must be disentangled. 

Recent theoretical papers show that breakup calculations are sensitive mainly to asymptotic properties of the 
projectile [ij, [l^, [l^. Both the ANC of the initial bound state, and the phase shifts in the continuum play significant 
roles in the dissociation. However, they do not prove the breakup not to be sensitive to the interior of the bound-state 
wave function. 

' In most reaction studies, peripherality of a reaction refers to insignificant contributions to the breakup cross section 
from small distances between the projectile and the target (e.g. [13, [3). However here we are interested in periph- 
erality relative to the internal coordinate of the projectile, the distance between the core and the valence nucleon. 
This notion of peripherality is intimately related to the sensitivity of the reaction to the interior of the projectile wave 
function. Since it is not directly related to the projectile-target peripherality, it deserves particular attention. 

To analyze the role played by the internal part of the wave function in the dissociation, we perform breakup 
calculations using two descriptions of the projectile that differ only in the interior of the wave function. To obtain 
such a pair of descriptions, we use the supersymmetric transformations developed by Baye p^ . [20| . Their great 
advantage is that the initial Hamiltonian and its supersymmetric partner exhibit identical asymptotic properties. 
Therefore, both descriptions exhibit the same phase shifts and ANC, but differ strongly in the interior. 
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The difference between calculations performed using both descriptions will tell us about the sensitivity to the interior 
of the projectile wave function. In this study, we consider two projectiles: a proton-halo nucleus (*B), and a neutron- 
halo nucleus (^^Be). The calculations are performed at intermediate (40-70 MeV/nucleon) and low (26 MeV) energies 
using either the dynamical eikonal model of breakup '^l', ^23\ or the Continuum Discretized Coupled Channel technique 
(CDCC) [l^ll^l. We also consider Coulomb and nuclear dominated reactions. In each case, the differences between 
both descriptions are analyzed through different observables (energy, angular, and parallel- momentum distributions). 

In Secini we present the projectile descriptions used in the calculations, including a summary on the supersymmetric 
transformations. The results of our analysis are displayed in Sec. IIIII The conclusions are drawn in the last section. 



II. PROJECTILE DESCRIPTIONS 

A. Supersymmetric transformations 

In both the dynamical eikonal [U and the CDCC [H, [l^l models of breakup reactions, the projectile is 
described as a two-body system: a spherical, structureless core c to which a pointlike fragment / is loosely bound. 
The internal structure of the projectile is modeled by the Hamiltonian 

Ho = -^A + V„ir), (1) 

where r is the relative coordinate of the fragment to the core, and fi is the reduced mass of the c-f system. The 
potential Vq describes the interaction between the core and the fragment. It is adjusted to reproduce the loosely 
bound states, and some of the resonances of the projectile. It contains both Coulomb and nuclear terms. For the 
present study, we choose Vq deep enough to obtain an additional deeply bound state in the partial wave of the physical 
ground state. This spurious state has no physical meaning, and merely serves to create a node in the interior of the 
wave function of the loosely bound state. 

To obtain a second description of the projectile, that differs from the first only in the interior, we remove the 
spurious deeply bound state. This eliminates the node in the wave function of the physical ground state. However, we 
want to preserve the ANC of the bound state as well as the phase shifts obtained with Vq, because they play a major 
role in breakup calculations [13, [H, [l^ . As shown by Baye [19l [20j , this can be done using a pair of supersymmetric 
transformations. The first removes the deeply bound state, without altering the other levels of the spectrum [2^ . 
The second restores the phase shifts of Vq that are modified by the first 1^, 2^ . The resulting potential V2 exhibits 



the same spectrum as the initial Vq but for the removed state. It is also phase equivalent as it gives the same phase 
shifts. Moreover, the ANC of the remaining bound states are identical in both descriptions. 

The supersymmetric elimination is performed only in the partial wave of the removed state. Therefore, V2 depends 
strongly on the angular momentum even if initially a single potential is chosen to describe all partial waves. In the Ij 
partial wave, the new potential reads [l^l 

Vi'ir) = V^'ir) - ^l^\nj\l{E,y)fdr', (2) 

where u^^ is the reduced radial wave function of the removed state of energy Eq. The potential V2 differs from Vq only 
in the range of Uy. Similarly, the wave functions of the Hamiltonian eigenstates are affected only at short distances. 
After the supersymmetric elimination, the wave functions read [20| 

where My and uf^ are the reduced radial wave functions obtained with Vq and V2, respectively, for either bound or 
scattering states. The asymptotic properties of Vq (ANC of the remaining bound states, and phase shifts in the 
continuum) are therefore preserved by these transformations. 



B. Descriptions of B and ^ Be 



Our description of ^B is based on the two-body model used in Refs. [2J,|2g|. The parameters of the ^Be-p potential 
are the same as in those former papers but for the depth of the central term in the p3/2 partial wave. For that 
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FIG. 1: Effective potentials describing the core-fragment interaction for *B (a), and ^^Be (b). Both the initial deep potentials 
(dotted Unes), and their supersymmetric partners (full lines) are depicted. Radial wave functions of the loosely bound state 
of *B (c), and ^^Be (d). Their asymptotic behavior extrapolated down to r = is also shown (dashed line). Radial integrand 
appearing in the first-order approximation of the transition from the ground state to a continuum state at _E = 1 MeV in the 
sl/2 partial wave for *B (e), and in the p3/2 partial wave for '^^Be (f) (see Sec. IIII B|l . 



partial wave we consider a much deeper Woods-Saxon form factor, which reproduces not only the loosely bound state 
of *B but also a spurious deeply bound state (at -Eop3/2 = —57.2 MeV). This potential (including the Coulomb and 
centrifugal terms), and the corresponding radial wave function of the physical ground state are displayed in Figs.[TJa) 
and (c), respectively (dotted lines). 

To obtain the second description of the nucleus, we remove the spurious p3/2 state using the supersymmetric 
transformations (see Sec. Ill A[) . This gives us the new potential in the p3/2 partial wave plotted as a full line in 
Fig. [Ija). The corresponding wave function of the physical p3/2 state is shown in Fig. [Ijc) (full line). As expected, 
both wave functions are identical above 4 fm, where they coincide with their asymptotic behavior (dashed line). 
However, they strongly differ in the interior, where one exhibits a node, while the other is nodeless. 

To describe ^^Be, we consider the same ^°Be-n potential as in Ref. [23|. This model includes a Pauli forbidden state 
in the sl/2 partial wave of the physical ground state (i?osi/2 = —32.7 MeV). To obtain the second ^^Be description, 
we simply remove that state using the supersymmetric transformations. The initial (dotted line) and supersymmetric 
(full line) potentials are displayed in Fig.[ljb). The wave functions of the physical loosely bound state are displayed 
in Fig.^d) for both descriptions. As in the ^B case, they differ at small distance, while they concur above 4 fm. 
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FIG. 2: Breakup cross section of *B on Pb at 44 MeV/nucleon as a function of the energy between the proton and the ^Be 
after dissociation. The results are obtained using either the deep potential (dotted line) or its supersymmetric partner (full 
line) . 



III. BREAKUP CALCULATION 



A. Coulomb breakup at intermediate energy 

Using the two descriptions of ^B, we perform calculations of its breakup on lead at 44 MeV/nucleon. It corresponds 
to the experiment performed at MSU by Davids et al. [28t . For this study, we use the dynamical eikonal approximation 
[2ll , ^[^ . To simulate the interactions between the projectile components and the lead target, we follow Mortimer et 

The results of the calculations are illustrated in Fig. [21 where the breakup cross section is plotted as a function of the 
relative energy between the ^Be core and the proton after dissociation. The difference between the two calculations is 
completely negligible (about 0.5 %). This confirms that breakup probes only the tail of the projectile wave function, 
since it is not sensitive to the large differences in the interior. 

A partial wave analysis of this cross section shows that the same result is obtained for the contribution of all partial 
waves. The tiny difference observed in Fig. [5] is not due to some accidental cancellation. Negligible differences are 
also observed in angular and parallel-momentum distributions. Note that the presence of the spurious deeply bound 
state in the p3/2 partial wave of the deep potential does not affect the breakup calculation. The cross section for the 
inelastic process leading to the population of this state is indeed of 8 /ib, which is five orders of magnitude below the 
total breakup cross section (0.874 b for the deep potential, and 0.869 b for its supersymmetric partner). 

The present result is in perfect agreement with the analysis presented in Ref. [29i|. In that paper, it was found 
that eliminating the forbidden bound states from the ^^Be spectrum has a very weak effect on the calculation of its 
breakup on lead at 70 MeV/nucleon. 



B. First-order analysis 

With the aim of understanding qualitatively this result, we make a first-order analysis. In this approximation [30| . 
the breakup is assumed to occur in one step from the initial bound state to the continuum. The reaction being 
Coulomb dominated, we consider here only the Coulomb interaction, which we expand into multipoles. The breakup 
cross section can then be expressed analytically as the sum of the contributions of the multipoles [s^]. The same 
result is obtained within this framework: the difference between both descriptions is negligible. 

The dependence on the projectile description of the contribution of the multipole E\ comes from the radial integral 

/>oo 

^eI^ ui,{E,ryui,,,{Eo,r)dr, (4) 

where uij is the radial wave function describing the c-f continuum in the Ij partial wave at energy E, and uigj„ is the 
radial wave function of the initial bound state. The integrand in Eq. ^ is displayed in Fig. \V[e) for the dominant 
El transition from the initial p3/2 bound state of ^B to the sl/2 continuum state computed a,t E = 1 MeV. This 
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FIG. 3: Breakup cross section of *B on Ni at 26 MeV as a function of the ^Be scattering angle after dissociation. The results 
are obtained using either the deep potential (dotted line) or its supersymmetric partner (full line). The difference is so small 
that both curves are superimposed. 

integrand for the El transition in '^-'^Be from the sl/2 bound state to the p3/2 continuum state at E = 1 MeV 
is depicted in Fig. [Ijf). For both nuclei, the major contribution to TZsij comes from radii above 4 fm, where the 
wave functions obtained with the deep and supersymmetric potentials coincide. The interior of the wave functions 
contributes very httle to the integral mainly because of the large tail of the wave function of the loosely bound state. 
Moreover, the vanishing behavior of the radial wave functions near the origin, multiplied by the r'^ factor, strongly 
reduces the contribution of the interior. This explains qualitatively why the internal part of the wave function has 
nearly no effect on breakup calculations. It also indicates that one cannot choose any form factor of the wave functions 
in the interior. For example, one could use the extrapolation down to r = of the asymptotic behavior of the wave 
functions [see dashed lines in Figs.[IJc) and (d)]. This leads to the dashed curves in Figs.[lje) and (f). In this case, 
the interior contribution is obviously overemphasized. If non-vanishing wave functions are used, it is necessary to 
resort to a cutoff at low radius, as Typel and Baur in Ref. [lB| . 



C. Coulomb breakup at low energy 



To see if the result of Sec. IIII Al holds at lower energy, we calculate the breakup of ^B on nickel at 26 MeV. It 
corresponds to the experiment performed at Notre Dame [3l|. This reaction has been analyzed within the CDCC 
framework in Ref. |2J] . For the present study, the calculations are performed in a similar way using the code FRESCO 

The angular distribution of the ^Be core after breakup is shown in Fig. [31 The difference between both ®B descrip- 
tions is even smaller than in the previous case: the two results are indistinguishable. Interestingly this is true for the 
entire angular range. In particular, even at large scattering angle, where the nuclear interaction dominates, no differ- 
ence can be observed between both ^B descriptions. The same conclusions are drawn from the energy distribution. 

This additional result strengthens the conclusion of Sec. IIII Al by extending its validity toward low energies. The 
dissociation of loosely bound nuclei therefore never seems to reflect the interior of the projectile. Only the tail of the 
wave function is probed by such reactions. 



D. Nuclear-dominated breakup 



Finally, we investigate the influence of the interior of the projectile wave function onto nuclear dominated reactions. 
We choose the breakup of ^^Be on carbon at 68 MeV/nucleon. This reaction has been measured at RIKEN by Fukuda 
et al. JjJ . It has been successfully described using the dynamical eikonal approximation ^22j, t27|i| , as well as the CDCC 
technique [s^]. We use the former model with the same inputs as in Ref. 

The breakup cross section is plotted in Fig. |4] as a function of the relative energy between the ^"Be core and 
the neutron after dissociation. The sensitivity to the potential choice is similar to the previous calculations. Both 
descriptions lead to nearly identical cross sections. The same conclusion can be drawn from angular and parallel- 
momentum distributions. This result shows that even in nuclear-dominated breakup, the calculations are not sensitive 
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FIG. 4: Breakup cross section of ^^Be on C at 68 MeV/nucleon as a function of the energy between the neutron and the ^''Be 
after dissociation. The results are obtained using either the deep potential (dotted line) or its supersymmetric partner (full 
line). 

to the internal part of the projectile wave function. As in the previous cases, only the tail of the wave function is 
probed. 

The only place where the two descriptions give slightly different results is in the peak at 1.3 MeV, where they differ 
by 8 %. The peak is due to the presence of a |^ resonance at that energy ^27^. Since the difference appears only in 
the c?5/2 partial wave, where the resonance is reproduced, we attribute this effect to the particular behavior of the 
continuum wave function. This wave function exhibits a huge peak near the origin due to the resonant nature of the 
state. In a first-order viewpoint, this peak increases significantly the contribution of the interior to the radial integral 
of Eq. ([4]). Consequently, a larger difference is observed. This effect is in agreement with the analysis of Ref. 29], 
where it has been observed that the excitation to the i bound state of ^^Be is more sensitive than the breakup to 
the removal of the Pauli forbidden states. In that case, both the initial and final states are bound. As in the present 
case, the low-r contribution is enhanced. Therefore, due to the shorter range of the final-state wave function, the 
excitation to bound or resonant states seems more sensitive to the interior of the wave function than the non-resonant 
breakup. 



As mentioned in the introduction, other theoretical studies refer to peripherality of reactions as insensitivity to 
small projectile-target distances [13, [3- For compact projectiles, no contributions from small impact parameters 
implies the insensitivity to the internal part of the projectile wave function. However, this may not be the case for 
loosely bound projectiles and reaction processes with multi-step mechanisms. 

Up to now, we have focused on peripherality relative to the internal coordinate of the projectile. In this section, 
we analyse the sensitivity of the breakup probability Pbu as a function of the impact parameter b. In semiclassical 
models, this breakup probability corresponds to the contribution of one trajectory to the total breakup cross section. 
This is illustrated in Fig. [51 

The upper part (a) corresponds to the breakup on lead at 44 MeV/nucleon for a relative energy E — 0.5 MeV 
between the ^Be core and the proton after dissociation. The breakup probability computed using the deep potential 
is displayed in dotted line. The full line corresponds to the relative difference between the results obtained with the 
deep potential and its supersymmetric partner. The first important point to notice is that the reaction occurs at very 
large impact parameters: most of the contributions occur beyond impact parameter of 15 fm. Secondly, in the region 
where there are contributions to the cross section, there is less than 1 % difference between the calculations using the 
deep projectile potential and the supersymmetric partner, confirming no sensitivity to the projectile wave function 
in the interior. This is a case where the reaction is peripheral relative to both, projectile-target and core-fragment 
distances. However, the two peripheralities seem disconnected. Beyond 15 fm, the relative difference between the 
deep and supersymmetric potentials is indeed very constant. It does not decrease with b: the sensitivity to the interior 
of the projectile description is the same at large projectile-target distances (6 ~ 100 fm) than at intermediate ones 
(6 - 15 fm). 



E. Comparison with projectile-target peripherality 
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FIG. 5: Relative difFerence between the deep and supersymmetric calculations in the breakup probability as a function of the 
impact parameter. The breakup probability obtained with the deep potential is shown as well, (a) Breakup of *B on Pb at 
44 MeV/nucleon for E = 0.5 MeV. (b) Breakup of "Be on C 68 MeV/nucleon for E = 0.5 MeV, and E = 1.274 MeV (i.e. 
l^-resonance energy). 



We perform the same analysis for the dissociation of ^^Be on carbon at 68 MeV/nucleon. The breakup probability, 
and relative difference between the two descriptions are displayed in Fig. [5]^b) for two relative energies E: 0.5 MeV, 
and 1.274 MeV, the energy of the |^ resonance. In this case, the dominant contribution to the breakup cross section 
comes from impact parameters around 6 = 5 fm. Expectedly, this nuclear dominated dissociation is less peripheral 
than Coulomb breakup. However, the difference between both descriptions remains negligible, except for the resonance 
energy, as seen in Sec. lIIID] Therefore, even though the reaction appears to be surface peaked in the projectile-target 
distance, it is still insensitive to the description of the projectile at small core-fragment distances. 



IV. CONCLUSION 



In this paper, the sensitivity of elastic breakup to the internal part of the projectile wave function is analyzed. We 
compare calculations of dissociation of loosely bound nuclei performed with two descriptions of the projectile, which 
differ only in the interior of the wave function. These descriptions are obtained from phase equivalent potentials 
constructed with supersymmetric transformations (lol . [20| . 

We start from an initial deep potential adjusted to reproduce, not only the loosely bound state of the projectile, 
but also a spurious deeply bound state. Due to the presence of that unphysical state, the projectile wave function 
exhibits a node near the origin. A second description is obtained by removing the spurious state. The wave function 
therefore no longer exhibits a node. This removal is achieved by supersymmetric transformations, which preserve the 
asymptotic properties of the initial potential [l^, (i.e. ANC and phase shifts). This enables us to analyze the 
sensitivity of the calculations only to the interior of the wave function. 

We study this effect in the dissociation of two halo nuclei: ^B, and ^^Be. Both Coulomb and nuclear induced 
reactions are considered, as well as intermediate and low energy regimes. In all the analyzed cases, no difference 
between the deep potential and its supersymmetric partner is obtained, although the corresponding wave functions 
of the projectile strongly differ in the interior (presence or absence of a node). This is observed in several observables 
such as energy, angular, and parallel-momentum distributions. We conclude that the elastic breakup of loosely bound 
nuclei is peripheral in the sense that it is not sensitive to the interior of the description of the projectile. Only the 
asymptotic properties of the projectile model are probed through these reactions; namely the ANC of the bound 
state, as suggested in Refs. and the phase shifts in the continuum Therefore, two descriptions 
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of the projectile corresponding to the same ANC but different spectroscopic factors will lead to the same theoretical 
prediction. This shows that spectroscopic factors extracted from breakup measurements will contain an error due 
to the imprecision of the projectile description at short distances, which is not probed by breakup reactions. The 
extraction of an ANC from breakup is less model dependent. Nevertheless, it will still contain uncertainties associated 
with the projectile-target interactions and the description of the projectile continuum. 
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